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We consider brane cosmology when the 4D Ricci scalar term is added to the 
5D Einstein-Hilbert action and discuss the role that the addition of this term has 
on the brane-bulk system. The induced brane dynamics is shown to be the usual 
I Einstein dynamics coupled to a modified energy-momentum tensor which is well 

defined once the 5D Einstein equations are solved in the bulk. The 5D Einstein 
equations valid everywhere in the bulk, but not in the brane, are projected on the 
brane. Then making use for the embedding of the brane in the bulk of the Israel 
junction conditions, modified by a source term coming from the addition of the 
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intrinsic curvature scalar in the bulk action, it is possible to obtain the effective 4D 
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Einstein equations on the brane consistent with the bulk geometry. 
PACS numbers: 04.20.Jb : 04.50.-h 



1. Introduction 

In the present work we investigate the cosmological evolution of the brane-bulk system in 
the framework of the Randall- Sundrum (AjdSs scenario [1]. The effective 4D gravitational 
equations in the brane without curvature correction terms were first obtained by Shiromizu, 
Maeda and Sasaki [2]. These equations have been later recovered and generalized both on 
the brane and in the bulk taking into account the effect of a general bulk energy-momentum 
tensor and either the asymmetric embedding [3] or the accelerations of normals [4] . However, 
even employing more generalized gravitational actions, the derived 4D Einstein equations 
do not in general form a closed system due to the presence of a Weyl term which can only 
be specified in terms of the bulk metric, so other equations are to be written down and 
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different procedures arise in splitting the non-Einsteinian terms between bulk and brane 
[5]. We assume omogeneity and isotropy in the three ordinary spatial dimensions but these 
symmetries cannot be extended to the extra dimension due to the presence of the brane, so 
all the physical quantities will depend on time and on the extra dimension. The solutions 
of the 5D Einstein equations that we shall obtain will be valid strictly in the bulk. Then we 
project on the brane at y = those equations and making use of the junction conditions, 
modified by an additional term coming from the curvature correction, as the boundary 
conditions imposed for the embedding of the brane in the bulk we obtain the effective 4D 
Einstein equations on the brane consistent with the bulk geometry. The method used here of 
deriving the 5D solution that is then projected onto the brane to study the brane dynamics 
was used extensively in the early days of braneworld cosmology by P.Kanti et al [6]. The 
paper is organised as follows. In the next Section, before giving a brief account of our 
method, we summarize the results obtained by Kofinas [7] when the term is included 
in the 5D action. Such a term, that was considered by a number of authors in the literature 
(see [8], [9], [10] and references therein), is generically introduced by quantum corrections 
coming from the bulk gravity and its coupling with matter confined to the brane, moreover 
its inclusion brings a convenient decomposition of the matter terms. In Section 3 we find 
the related equations in the brane, assumed infinitely thin and Z 2 symmetric in the bulk. 
In Section 4 we use the flexibility of the 5D solution to describe some cosmological models 
in 4D. Finally, in the Appendix we show how the 5D dynamical solution we adopt can be 
obtained starting from a 5D static solution. 

Conventions. Throughout the paper the 5D metric signature is taken to be (+, +, +, —,e) 
where e can be +1 or —1 depending on whether the extra dimension is spacelike or timelike, 
while the choice of the 4D metric signature is (+, +, +, — ). The spacetimes coordinates are 
labelled x l = (r,$,ip),x 4 = t. The extra coordinate is x 5 = y. Bulk indices will be denoted 
by capital Latin letters and brane indices by lower Greek letters. In what follows we choose 
units such that H — c — 1. 

2. Braneworld Einstein field equations 

In this section we recall the results obtained by Kofinas [7], which we shall use in the 
following, giving a brief account of their derivation. Once we have solved the equations in 
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the bulk, the form of the induced equations will allow us finding brane solutions following the 
methods of General Relativity with a well defined energy- momentum tensor. The starting 
point in [7] is a three-dimensional brane E embedded in a five- dimensional spacetime M. 
For convenience the coordinate y is chosen such that the hypersurface y = coincides with 
the brane. The total action for the system is taken to be 

S =7T2 f \f Z ^g( {5) R-2A 5 )d 5 x + -^ [ v /r ^( (4) ^-2A 4 )rf 4 x 

^ K 5 JM * K A JY, 

+ / v 7 -^ ^gL™ at d 5 x+ [ v 7 - WgL™ at d A x (1) 

JM JY 

The constants k\ and k\ are given by 

K l = 8nG 5 = M 5 " 3 , k\ = 87rG 4 = Mf (2) 

where M 5 and M 4 are the Planck masses. Varying (1) with respect to the bulk metric Qab 
one obtains the equations 

®G B A = -A 5 5 B + 4 (^T B T B 5(y)) (3) 

where 

{loc) n = -~ 2 J^f ( i4) G B A - «2 (4) T? + A 4 ^) (4) 

is the localized energy-momentum tensor of the brane. ^Gab and ^Gab denote the Ein- 
stein tensors constructed from the bulk and the brane metrics respectively, while the tensor 
h-AB = 9ab — suaUb is the induced metric on the hypersurfaces y = constant, with n A the 
normal unit vector on these 

<-A 

n A = f, n A = (0,0,0,0,5$) (5) 

The scalar $ which normalizes n A is known as the lapse function and in a cosmological 
scenario which we shall consider later, it will depend on t and y. The way the coordinate y 
has been chosen allows to write the five- dimensional line element, at least in the neighborood 
of the brane, as 

dS 2 = ()ab dx A dx B = dx»dx v + e <5> 2 dy 2 (6) 



4 



Using the methods of canonical analysis [11] the Einstein eqs. (3) in the bulk are split into 
the following sets of equations 

K;., v -K,, =e4^) T y (7a) 
K»K» - K 2 + e (4) i? = 2 e (A 5 - 4 ^T%) (7b) 
v - + $ K K» - e $ (4) i# + £ A = MM T ^ - - {loc) T SS\ 5(y) 



- e k\ $ (5) T^ + - $ ( K 2 ( 5 )t - 2A 5 ) ^ (7c) 



<9y 

where fT^ is the extrinsic curvature of the hypersurfaces y = constant: 

K " = J*W' Ka ' = (8) 

The Israel junction conditions [12] for the singular part in eq. (7c) are 

[K^] = -e 4 $ (V) T A. _ 1 (ioc) T 5 fj (g) 

where the square brackets mean discontinuity of the quantity across y = and $ represents 
<3> at y = 0. Consequently, considering a Z 2 symmetry on reflection around the brane, (9) 
becomes 

(4) G: = -A 4 K + ^ (4) T^ + - (^ - 7f<f£) (10) 

where K^ v = K^(y = + ) = — K£(y = 0~) and r c = 4/4 * s a crossover term which 
determines the region of validity of conventional four-dimensional General Relativity. From 
eq. (7a) it follows that the tensor satisfies the conservation law ^Tff = provided 
(5)Tw = 0, which means no exchange of energy between brane and bulk. The quantities 
are still undetermined and should be obtained from some exact solution of the global 
five-dimensional spacetime. To determine the equations on the brane one can follow the 
method suggested in ref. [2], but this will reveal a difficult task due to the necessity of 
taking into account the evolution of the Weyl term to close the system of equations. A 
different approach, as discussed by Binetruy et al. [13], is to solve the 5D Einstein equations 
strictly in the bulk (y ^ 0) and then to take the brane into account by using the Israel 
junction conditions. In this work we shall keep in mind this latter approach but, having 
added the term, we shall make a different use of the junction conditions. More in detail, 
we start projecting on the brane the solution obtained in the bulk without considering the 
distributional part at y — 0. This will be done using the geometrical identity 

iA) R A BCD = (5) Rnkl h A M h N B h«h L D + e [K* K BD - K B K BC ) (11) 
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and taking suitable contractions from the above relation. So it is possible to construct 
the four and five- dimensional Einstein tensors and to get finally the parallel to the brane 
equations 

+^ - *J + | - **) K - <T ^ (12) 

Here is the "electric" part of the bulk Weyl tensor, while T and C are the limiting 

values of those quantities at y = + or 0~. Once we have solved the Einstein equations 
strictly in the bulk we can make explicit the various terms appearing in the right-hand side 
of (12). Now we impose to the above solution boundary conditions to take into account 
the physical presence of the brane. This can be done if we consider equations (10) as the 
boundary conditions imposed for the embedding of the brane in the bulk so we have to 
equate the two independent equations (10) and (12). In this way, however, we would obtain 
the "bare" quantities A 4 and ^Ttf but not the effective quantities as seen by an observer 
confined to the brane. For a brane observer eqs. (10) are instead written as the usual 
Einstein equations 

^G» = -A 4 , eff 5i: + Kl ( (4) T^) e// (13) 

So the cosmological constant A 4je // and the effective energy-momentum tensors {^Tjf) e ^ 
can be obtained equating the right-hand sides of eqs. (12) and (13). 

3. Dynamics in the brane-bulk system 

We consider the 5D metric in the form commonly used in cosmological applications 

dS 2 = a 2 (t, y) da\ - n 2 (t, y) dt 2 + e $ 2 (t, y) dy 2 (14) 

where 

da l = J^— + r \ d tf + sin 2 (is) 

and k = +1,0,-1 is the curvature index. Having specified the form of the metric, we 
now turn to the 5D Einstein equations (3) considered strictly in the bulk, that is, without 
the energy- momentum tensor at y — 0. These equations can be solved once the structure 
and the content of the bulk come as a result of a physically acceptable theory in higher 
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dimensions. Exact time-dependent solutions which generalize the static solutions were con- 
structed using diffeomorphism invariance [14,15]. Kehagias and Tamvakis [16] transformed 
the static Randall-Sundrum (A)dS$ solution into a dynamical one by considering boosts 
along the fifth dimension and found a time-dependent 5D solution for a bulk with vacuum 
energy but otherwise empty and with vanishing Weyl tensor. In the present work we want 
instead to consider (see (12)) a bulk where the non- localized energy-momentum tensor and 
the "electric" part of the Weyl tensor are different from zero, so we should start by a well 
defined bulk matter content described by the tensor and then solve the field equations. 
However our aim is to overcome the problem of the brane field equations being non-closed 
so, to give an illustrative example of our method, we shall proceed in a bit unhorthodox 
way. In order to have a simple and non-trivial dynamical 5D solution we start from a static 
Randall-Sundrum (A)dS 5 bulk and we construct, generalizing the transformations in [16] 
a dynamical 5D line element with non-vanishing 5D Weyl tensor. Subsequently we obtain 
the correspondent energy-momentum content using the Einstein equations. This manner 
of proceeding may be justified by the fact that this work is mainly focused on the brane 
phenomenology of the model. Here we anticipate the main results of the procedure and de- 
fer to the Appendix for detailed calculations. Our dynamical line element will be obtained 
transforming the static Randall-Sundrum (A)dS 5 metric, where the three-space is not nec- 
essarily flat but has a curvature index k = +1,0, —1, into a dynamical one by considering 
boosts along the fifth dimension. Then we take into account the Einstein equations in the 
bulk, away from the brane at y — 0, and require that there is no energy flow from the brane 
towards the bulk and vice-versa, which implies = 0. The above constraint is easily 

satisfied if one chooses wave-like expression for the metric coefficients so, assuming the Z 2 
symmetry y — > — y, it follows that a(t, y), n(t, y) and $(t, y) in the line element (14) will be 
function of w = t — \ \y\ with A a dimensionless constant. Metric coefficients in the form of 
plane waves propagating in the fifth dimension have previously been used in the literature 
in somewhat different contexts [17,18,19,20]. Finally, the bulk line element away from the 
brane was found to depend only by the scale factor a(t, y) = a(t — \ \y\) in the form 



dS 2 = a 2 da\ — 



2 7 2 A 2 



1 



( 



K 2 a 2 + y / /t 4 a 4 + 4£ 7 2 A 2 (a) 2 ) dt 2 




(16) 
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where k is the constant scale factor for the extra dimension, the superscribed asterisk 
denotes derivative with respect to w and 7 is a dimensionless constant which comes from the 
constraint = 0, namely a = 711 $. It should be noted that a(t — A \y\) = corresponds 
to a scale factor singularity for the 5D model which is similar to those that occur in the 4D 
Friedmann-Robertson- Walker models. From the following curvature invariants 

®R=-20e K *+ 2(f 6 \, ^R AB ^R AB = m^- + ™* , (17) 

a\t-Mv\) a 2 (t-X\y\) a A (t-\\y\) 

^Rabcd^R ABCD = 40« 4 - 2 l 4€k f, + 12fc2 (18) 

a 2 (* — A |y|) a 4 (t-A|y|) 

we see that there is no other singularity except the one which may unaivodably occur if the 
scale factor of the fifth dimension vanishes in some (t, y) hyperplane. 

Now we calculate the 5D Einstein tensor and from the field equations considered strictly in 
the bulk we obtain the cosmological constant A 5 and the energy-momentum tensor as 

A 5 = -Qek 2 , <®T% = diag(p B ,PB,PB,-PB,P±) (19) 

where 

4pb = ~ 2(+ ~ I rr , *4 p B = -K 2 5 p ± = 3 ^ (20) 
a(t-\\y\) a 2 (t-X\y\) 

the subscript b referring to the bulk. A comment is needed about the cosmological fluid 
described by the energy-momentum tensor which arises from the curvature index k in eq. 
(15). It obeys an equation of state p B = (7 — 1) p B with barotropic index 7 = 2/3, its 
pressure and energy density are proportional to k and scale as a~ 2 . All these features lead 
to the interesting possibility that the energy-tensor (20) can describe a fluid composed of 
cosmic strings, as discussed by a number of authors [21]. We can use the flexibility of the 
metric (16) to choose many different 5D scale factors, but clearly each choice must meet the 
necessary requirements to give models acceptable on physical grounds. 
Let us now deal with the brane dynamics. We consider homogeneous and isotropic geome- 
tries in the brane so the effective tensor {^Tff) e ^ will describe a cosmological fluid endowed 
with pressure p e ff and energy density p e //. Equating eqs. (12) and (13) we obtain 

2 A 5 e (\ da\ {Ida 2 dn\ k . . 

-A 4 ,„ + Kl p,„ = - T - 5 ^--j^-- + --j r - I (21) 

a 2 A 5 3e (\ da\ 3k . . 

- A 4 , e // - «4 Peff = —2-& {adi) --% (22) 
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The effective cosmological constant A 4i6 jj may be different from A 5 /2 = — 3e/t 2 , its value 
being modified by possible additive constant terms contained in (21) and (22). The Einstein 
tensor which appears in the left-hand sides of (12) and (13) is constructed from the 

brane metrics 

ds 2 = a 2 (t)da 2 k -n 2 (t)dt 2 (23) 

Now, in higher-dimensional theories there is the question of which metric frame is the correct 
representation of our four- dimensional spacetime. In many braneworld theories, the physical 
metric in 4D is identified with the induced one, while in other approaches the physical 
metric either is assumed to be conformally related to the induced one or is determined by 
the condition of classical confinement in the absence of non-gravitational forces [22]. For 
the sake of simplicity, here we choose to identify the metric (23) with the induced one, 
so we have a(t) = a(t,0) = a (t) and n(t) = n(t,0) = n (t). It follows that eqs. (12) 
and (13) are identically satisfied by the Einstein tensor ^G^ constructed from (23). From 
the knowledge ao(t) and no(t) one can also obtain other cosmological quantities such as 
the Hubble parameter H = a'o/(noao) or the deceleration parameter q = — (aod'o)/Gk) 2 + 
(do rio) I (do no) . Difficulties may instead arise from the exact evaluation of the 4D proper 
time r when dealing with the integral r — Jn dt and a generic value of n (i). Finally, if 



we define 



k 3k 

k\P4> = K>lp e ff + ~2 and K l p<t, = k\ p e ff j (24) 



a a 



we can model the fluid in terms of a scalar field 0, minimally coupled to Einstein gravity 
and self- interacting through a potential V(<f>), with pressure and energy density given by 



V4> = ±<t>o/{2nl)-V (25) 



• 2 



p = ±0o /(2n 2 ) + V (26) 
where the upper (lower) sign corresponds to a standard (phantom) scalar field. 



4. Some possible brane scenarios 

In this section we describe two of the possible brane scenarios consistent with our bulk 
solution. We shall choose simple values for the scale factor a(t, y) and then determine A^// 
and (^T^) together with the parameters q and H. In the following we shall consider a 
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spacelike fifth dimension, so e = 1 and give relevant 4D quantities as a fuction of the 4D 
proper time r. 

A) Let us first consider the case a(t,y) = ( , jX/k) sin/t(t — X\y\). 

This choice gives no(t) = 1 so the coordinate time t now coincides with the 4D proper time 
r and therefore 

a o( r ) = (iX/n) sin k r, q(r) = tan 2 k t , H(t) — k cot k t (27) 

The evolution of the universe begins with a big bang at r = 0, reaches a maximum (a ) max = 
(7 A)/k and terminates with a big rip at kt — ir. We do not give numerical values for 7 
and A while k depends on the scale factor of the fifth dimension. The cosmological constant 
A 4 , e // and the tensor ( (4) T^) ( are given by 

a „ 2 (7 2 A 2 + k) 3 ( 7 2 A 2 + k) 

A 4,e// = - 3 K , p e // = - 9 2/^ ' P*ff = K 2 n 2(^ ( 2S > 

K 4 a {T) K 4 a yr) 

As to the standard scalar field, we have: 

0= — In tan — , V = — = , V U ) = cosh 2 29 

/t 4 2 re| sin 2 kt «| ^2 

The same results, starting from different points of view, were obtained in [23]. 

B) Now let us consider the case a(t,y) = (7A//C) (1 — n(t — A 1 2/ 1 ) 1 . 



This choice gives n (t) = \/2/(y/5— 1) (1 — Kt) 1 so the relation between the coordinate 



time t and the 4D proper time r is (1 — Kt) — exp 



l)/2 kt 



Therefore 



Qo(t) = — exp ( \l n 1 kt \, q = -l, H = \l v " - « (30) 



ft 




The cosmological constant A 4 e jj and the tensor {^T^ e ^ are given by 

3 (>/5 — 1) 2 k 3k <o-x\ 

A4 ' e// = 2 Pe// -~^I(7)' Pe// -^(7) (31) 

The values of a(t, ?/) chosen in the previous illustrative examples reproduce results already 
known in the literature. Less simple choices for the 5D scale factor may describe new 
brane scenarios but also may require a more involved treatment. In conclusion, this paper 
investigates the influence of the term included in the bulk action on the spherically 
symmetric braneworld solutions. The brane dynamics is made closed by using the modified 
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junction conditions as the boundary conditions for the embedding of the brane in the bulk, 
so it is possible to obtain brane cosmological solutions consistent with the bulk geometry. 
We started from a particularly simple time-dependent solution in the bulk away from the 
brane, but other physically acceptable solutions in the bulk can be considered provided that 
the related brane dynamics is in accordance to the observations on the brane. 



APPENDIX: Transforming static bulk solutions into dynamical ones 



The Randall-Sundrum (AjdSs model is the simple braneworld with curved extra dimen- 
sion that allows for a meaningfull approach to cosmology, therefore we start from this model 
but, at this point, we do not yet require Z 2 symmetry on reflection around the value Y = 
so we write 

dS 2 = e -2nY ( A l da l- dT l) +£dY 2 (A.l) 

Here k and A are, respectively, the constant scale factors for the extra dimension Y and for 
the ordinary three-space and da\ is the line element of maximally symmetric three-spaces 
with curvature index k = +1, 0, —1: 

dr 2 



^ = 1 - kr- 



+ r 2 (A? 2 + sin 2 tidip 2 ) 



(A.2) 



Since our purpose is to describe the time evolution on the braneworld, we need to transform 
the static bulk solution (A.l) into a dynamical one. This goal was already achieved in 
literature where dynamical solutions are derived from the static Randall-Sundrum (A)dSs 
metric by considering boosts along the fifth dimension [16]. Applied to the actual case, we 
generalize those transformations as 

X 



X 



K 



(A.3) 



kY 



F(t,y) + G(t,y) 



X 



K 



where F(t, y) and G(t, y) are dimensionless functions and x is a constant, with the dimensions 
of k, related to the boost along the fifth dimension. The coordinate y is chosen so that the 
hypersurface y = coincides with the brane. 
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As a result the metric (A.l) becomes: 
1 



dS z 



(F + G-iy 



« -x 



A 2 da\ + 



k 2 -ex 2 
k±x 2 



[- k 2 (dF) 2 + ex 2 {dG) 2 } 



(A.4) 

Note that the static line element (A.l) can be recovered from the above equations on con- 
dition that as x —> it results F & 1 — xt — e (x 2 /k 2 ) e Ky and G « xt + e Ky . 
The line element (A.4) is in the form commonly used in cosmological applications 

dS 2 = a 2 (t, y) da\ - n 2 (t, y) dt 2 + e $ 2 (t, y) dy 2 (A.5) 

Now we can choose suitable functions F and G to obtain explicit expressions for the metric 
coefficients a, n and <£>. Comparing eqs. (A.4) and (A.5) we get 

A 2 {k 2 - x 2 ) 



a 2 (k 2 - e x 2 ) 
A 2 k 2 x 2 {k 2 — x 2 ) 

a 2 (K 2 -ex 2 ) 
A 2 k 2 x 2 (« 2 - X 2 ) 



K 



K 



dF 
~dt 

dF 
dy 



k 2 (F + G-l) 2 

2 / ^v\ 2 



- £X 



dG 
~dt 

dG 
dy 



n 



K 2 I 



fdF_\ (dF_ 
\dt) \dy~ 



dt J \dy 



(A.6a) 
(A.6b) 

(A.6c) 
(A.6d) 



Once the new metric coefficients are known it is possible from the Einstein equations (3) 
in the bulk, that is, away from the brane at y — 0, to obtain the energy-momentum tensor 
(5)j>b_ rj^jg can achieved by recalling that in the coordinate system (A.5) the non- 
vanishing components of the Einstein tensor G^ are 



G v 



G y 

y 



rr 



e 

+ $2 



n z 



$ Id $ /2d 
$ a $ \ a 



+ 



a a 



a \a 



2n\ 
n J 




a' 



— + - - + 



a: 



2n' 



n 



+ 



3e fa" 

$2 



.'2 



+ 



$' /2a' n' 
$ \ a n 

a'$' 



a 



a^ 

,'2 



a a" dn\ 3e /a' 
s a a 2 any $ 2 \ a 2 



a'n' 



an 



a 

3A; 
'a 1 



3k 



3e 

$2 



o c i a an u >f 

an 




(A.7a) 
(A.7b) 

(A.7c) 
(A.7d) 



Here a dot and a prime denote partial derivatives with respect to t and y, respectively. In this 
work we require that there is no energy flow from the brane towards the bulk and vice- versa 
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so it must be ( 5 )Tf = 0, therefore the choice of the functions F and G must give accordingly 
^G y t = 0. However, as eq. (A.7d) shows, there is no energy flow only for suitable values 
of the metric coefficients. A particularly simple choice which makes ^G v t = is to assume 
that the metric coefficients in the bulk have the form of plane waves propagating in the fifth 
dimension, so they become functions either of the argument u = t — A y or of the argument 
v — t + A y. Of course the particular metric which we finally obtain is dependent on this 
choice. In detail, from ^G\ = we can derive 

1 da(u) 1 da(v) 



7 



(A.8) 



n(u) du n(v) $(i>) dv 

where 7 is a dimensionless constant. Now we shall assume the Z 2 symmetry y — > — y and 
construct a solution of eqs. (A. 6) by matching a solution depending only on u (for y > 0) 
to a solution depending only on v (for y < 0). 

The result in (A.8) suggests to multiply (A. 6b) by (A. 6c) so, taking into account (A.6d), we 
have 



a 2 (K 2 -ex 2 ) 



dF\ (dG_ 
~dt) \dy~ 



dF\ (dG 



dy J \dt 



A 2 ( K 2 _ X 2) 

Eliminating G by (A. 6a) and taking the square root one finally obtains 

da 



n 2 & 



k 2 -ex 2 



dF 
~dt 



, < + (if) (I) 

dyj \oy J \dtj 



n $ 



(A.9) 



(A.10) 



Ak 2 x\/k 2 - x 2 

Let us first begin working on the y > side. We put a(t,y) = a(t — \ y) into (A.10) and 
recalling (A.8) we obtain the following partial differential equation for F 



dF dF _ Ak 2 xV k2 - X 2 
dt dy 7 (k 2 — e x 2 ) 



(A.ll) 



The general solution is 



F(t,y) 



Ak 2 x\/ k 2 - X 2 



(* + Ay) + / ( _)(* -Ay) 



(A.12) 



2 7 A (K 2 -e X 2 ) 

The function /(_)(w) can be determined by (A.6d) after eliminating G by (A. 6a). The result 



is 



du a z (K z — ex ) K \du J 7 A y \du J 



(A.13) 



which can be integrated, once the scale factor a has been fixed, reminding that in the limit 
X — > it must be F(t,y) — > 1 and so also f(-)(u) — > 1. Of course if one is only interested 
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in determining n and $ from (A. 6b) and (A. 6c) it is sufficient the simple knowledge of the 
derivative of /(_)(«). Obviously n and $ will be a function of the scale factor a. Proceeding 
in an analogous manner when working on the y < side, we obtain 



where fr + )(v) and /(_)(«) are the same function / of the two different arguments v and u. 
As a consequence we can write the function F(t, y) on both sides of the brane at y = 
simply as 

F(t, y) = ~ % (t + A \y\) + f(t - A \y\) (A.15) 

27A (/t 2 - £X ) 

The function G(t,y) can then be easily derived from eq.(A.6a). Finally, we can obtain from 
eqs. (A. 6b) and (A. 6c) the metric coefficients n(t — A \y\) and $(t — A \y\) which are given 
as a function of a(t — A \y\) by 



n 2 = _L_ ( K V + V^ 4 a 4 + 4e 7 2 A 2 (a) 2 ) (A.16) 



$ 2 = __ k u + y « 4 a 4 + 4 e 7 2 A 2 (a) 2 J (A. 17) 

where the superscribed asterisk * denote derivative with respect to w — (t — A \y\) The bulk 
line element away from the brane is therefore 



dS 2 = a 2 da\ - ( K 2 a 2 + V K 4 a 4 + 4 e 7 2 A 2 (a) 2 ) df 

1 



kV + V^ 4 « 4 + 4e 7 2 A 2 (a) 2 V (A.18) 



2 7 2 

as given previously in eq. (16). From eqs. (A. 7) we get 

CT r = Gt = G% = 6sk 2 - \ (A.19a) 

Ik 

G\ = Gl= Qe K, 2 — — — — (A. 19b) 

in accordance with the 5D cosmological constant and the 5D energy-momentum tensor given 
previously in eqs. (19) and (20). 
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